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The two-electron orbitals for the U center have been computed numerically from the Hartree- 
Fock-Slater (HFS) equations in the point-ion lattice potential. The lattice relaxation of the nearest- 
neighbor ions is included in the model. The five lowest-lying cV-center states for NaCl, KC1, CdF2, CaF2, 
SrF-2, and BaF2 are given. The low-lying singlet states have the following order for increasing values 
of the energy: 'SCls, Is), ^(ls, 2p), and ^(ls, 2p) . The energy levels for the triplet states 3 S(ls, 2s), 
and 3 P(ls, 2p) lie between the energy levels for the ^(ls, Is) and 'PUs, 2p) states. The ordering of 
the triplet states depends upon the host crystal and the lattice relaxation. The predictions based upon 
the numerical HFS wave functions are compared with the predictions based upon past variational wave 
functions and with experiment. 
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1. Introduction 

The U center in ionic crystals consists of a nega- 
tive-hydrogen ion H", located substitutionally at an 
anion site. Because the H~ ion has a different mass, 
polarizability, and interaction with its neighboring 
ions than has the anion which it replaces, the V center 
changes both the phonon spectrum and the optical 
properties of its host crystal. Changes in the phonon 
spectrum include new features at energies substan- 
tially greater than the energies of the optic phonons 
of the host lattice [1-11]. 1 Changes in the optical 
properties include a characteristic electronic absorp- 
tion band, the U band, in the ultraviolet region [12-20] 
and are the concerns of this paper. 

The author continues in this paper his calculations, 
reported in part in [21], on the several states of the U 
center in the two alkali halides, NaCl and KC1, and in 
the four alkaline earth fluorides, CdF2, CaF-2, SrF2, 
and BaF2. Reference [21] contains discussions which 
stress the U band transition energies between the sin- 
glet ground and first excited states. The author now 
reports his results on the spatial extent of the ^-center 



1 Figures in brackets indicate the literature references at the end of this paper. 



orbitals, on the transition energies between the first 
and second excited singlet states, and on the transition 
energies between the two lowest-lying truplet states. 
A point-ion lattice model is used which incorporates the 
Hartree-Fock-Slater (HFS) procedure to compute 
numerically the defect electron orbitals. The model 
includes estimates for the correlation energy of the 
singlet states [22], for the exchange energy of the trip- 
let states [23], and for the lattice relaxation of the 
nearest-neighbor ions. The model also obeys com- 
pletely the Franck-Condon principle that the lattice 
relaxation does not respond to rapid changes in the 
electronic state of the U center when the U center 
undergoes an optical transition. This means that all 
low-lying f/-center states should have spatially compact 
electron orbitals, if the predictions of the model are 
to be internally consistent with its assumptions. 

The present point-ion model has some limitations. 
Its numerical wave functions are not orthogonal to 
the electronic states of neighboring ions and thereby 
do not account properly for the finite size of the neigh- 
boring ions, particularly when the lattice relaxes. An 
extended-ion model which uses numerical wave 
functions instead of variational wave functions for 
the defect electron orbitals has not been reported. 
Such a calculation would require an amount of com- 
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puter time which is substantially greater than that 
required by Wood and Opik [18]. These authors solved 
extended-ion models by using trial wave functions with 
two variational parameters and by devising accurate 
approximations for efficiently computing Coulomb and 
exchange integrals. Their methods use Slater-type 
functions and are not applied readily to the present 
numerical procedures. 

Keeping the above limitations in mind, we shall 
calculate here within the framework of the above 
model the energy levels of the three lowest-lying 
singlet states and the two lowest-lying triplet states. 
In some cases, these triplet states are degenerate or 
very nearly degenerate. 

Because expressions for the change in energy of a 
classical point-ion lattice with lattice relaxation of the 
neighboring ions and expressions for the total energy 
of the U center when the two electrons are in singlet 
and triplet states appear in sections 2 and 3 of [21], 
they are not repeated here. In section 2 of the present 
paper, the equations for absorption are given. Section 
3 contains a discussion of the results and some remarks 
on the validity of point-ion models. Finally, the ap- 
pendix contains a summary of the many assumptions 
upon which the HFS equations are based. 



r n i[s,ri(nl;n , r);<j] = rr s JR„i(r)r s R nl (r)47Tr 2 dr 

and 

r n 'v[s; y(nl; n'V)\ cr] 

= rT*$R n >v{r)T*Rnn>{r)ton*dr, 



(2) 



(3) 



where 5 = 1 or 5 = 3 and n is the nearest-neighbor 
distance of the perfect lattice. The principal and orbital 
quantum numbers for the radial wave functions R n i (r) 
are respectively n and /. The ratio, 



r e (nl) = [r n / (5 = 3)/r„/ (5=1)], 



(4) 



also indicates to what extent the radial functions have 
extended tails. The values of r e < 1 indicate compact 
states, and values of r e > 1 indicate diffuse states. 

3. Results and Conclusions 

In this section, the results of the point-ion model 
with lattice relaxation are reported. Table 1 contains 
the values of the input data which have been used. 



2. Absorption 

Because no emission bands have been associated 
with the U center [18], we shall study only the optical 
absorption which the defect electrons may undergo. 

The initial state of an optical transition is a relaxed 
state, 1 77;; cr, >. The symbol n denotes an electronic 
configuration of the two defect electrons which 
transforms as one of the following states, ^(ls, I5), 
3 S(l5, 25), 3 P(l5, 2p), ^(15, 2p), and ^(ls, 2s). The 
quantity cr characterizes the lattice relaxation; i.e., 
the n nearest-neighbor ions of the U center may move 
radially to the sites ~v' v = r v (1 — cr) for 1 ^ v ^ n. All 
other ions remain at their perfect lattice sites r„. We 
minimize E T (rji; cr) with respect to the lattice relaxation 
cr to obtain the energy of the initial state E x — E T (rji; cr,-). 
The value cr, is that value of cr for which E T (r)i; cr,-) 
attains its minimum value. The U center then under- 
goes an optical transition to an unrelaxed state \rjf; 
cr, >, which has an electronic wave function calculated 
for the same distortion cr, as that for the initial state 
|i7i; cr, >. The total energy of the final state is E/=E T 
(if, °"»)- The optical absorption energy for a transition 
from state i to state /is 



E(i,f) = E T (r)/, o-,) — E T (vt; cr/), 



(1) 



where each term on the right-hand side of eq (1) is a 
negative number. 

The expectation value of a given power of the radial 
coordinate r gives us information on the spatial extent 
of the two-defect-electron wave functions. As a measure 
of the spatial extent, the author chooses to consider the 
first and third powers of r for each of the orbitals; 
namely, 



Table 1. Input data for the point-ion model of the U center with 
ionic polarization, exchange energy, and Coulomb-correlation 
energy. The Pauling factor of the ^th and /xth ions is j3„ M . The ionic 
radius of the cation is p+ and the anion is P-. The quantity p, 
having the units of length, is the stiffness factor in the empirical 
Born-Mayer exponential form which characterizes the repulsive 
energy between the ^th and /^th ions. The Madelung potential 
constant at the anion site is a.v/. The quantity f\ is the nearest- 
neighbor distance (anion-cation) for the NaCl structure and is 
the lattice constant (cation-cation) for the CaF2 structure. The 
series coefficients C4, C«, andCs appear in the expansion in powers 
of the lattice distortion cr for the change in electrostatic energy E\ 
which occurs when a cation moves in the background of a perfect 
point-ion lattice potential; namely, E\——{n\f\) (C 4 cr 4 + c7 6 a- 6 
+ C«cr 8 -(- . . .), where n is the number of nearest-neighbor 
cations to an anion. The quantities j3 + + ' /3+-' and j3--' and olm 
are dimensionless. All other quantities are expressed in terms of 
atomic units (1 a.u. of energy = 27.2 eV = 4.36 attojoules and 1 a.u. 
of length = 0.0529 nm). 





NaCl 


KC1 


CdF 2 


CaF 2 


SrF 2 


BaF 2 


0++ 


1.25 


1.25 


1.50 


1.50 


1.50 


1.50 


0+- 


1.00 


1.00 


1.125 


1.125 


1.125 


1.125 


0- 


0.75 


0.75 


0.75 


0.75 


0.75 


0.75 


P 


a 0.599 


a 0.637 


0.582 


b 0.546 


b 0.560 


b 0.582 


P+ 


a 2.21 


a 2.77 


2.76 


c 2.21 


c 2.48 


c 2.76 


P- 


a 3.00 


a 3.00 


1.98 


c 1.98 


c 1.98 


c 1.98 


ocm 


1.748 


1.748 


4.071 


4.071 


4.071 


4.071 


f\ 


a 5.31 


a 5.93 


10.21 


c 10.32 


c 10.95 


c 11.71 


c 4 


d 3.579 


d 3.579 


e 1.865 


e 1.865 


e 1.865 


e 1.865 


c« 


d 0.9895 


d 0.9895 










c 8 


d 2.942 


d 2.942 











a M. P. Tosi, in Solid State Physics, edited by F. Seitz and D. 
Turnbull (Academic Press Inc., New York, 1964), Vol. XVI, p. 52. 

b A. D. Franklin (private communication). 

c G. C. Benson and E. Dempsey, Proc. Roy. Soc. (London) A266, 
344 (1962). 

d A. Scholz, Phys, Status Solidi 7, 973 (1964). 

e H. S. Bennett, J. Res. Nat. Bur. Stand. (U.S.), 72A (Phys. and 
Chem.), No. 5, 471-474 (Sept -Oct. 1968). 
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The numerical procedures give the energy eigen- 
values to an accuracy of \&E/E\ —0.001. They also 
give the self-consistent potential which appears in the 
HFS equations to an accuracy of | A V\V | ~ 0.01 , within 
the context of the point-ion model. The quantities 
A£ and AF are, respectively, the changes in the 
trial eigenvalue and self-consistent potential which 
occur between two successive iterations in the nu- 
merical itegration procedure. The numerical details 
on the spatial extent of the five lowest-lying states of 
the U center are reported only for KC1 and CaF2 with 
no lattice relaxation. These two crystals are quan- 
titatively representative of the remaining four crystals 
examined in this paper and illustrate respectively most 
of the salient features common to the NaCl-type and 
CaF2-type structures. Succeeding tables do contain, 
however, the energy levels for the five lowest-lying 
states in all six crystals. The computed data for the 
case in which the lattice relaxation and the electron 
orbitals are determined in a self-consistent manner are 
to within about 30 percent qualitatively similar to the 
computed data presented in tables 2 and 3. The only 
exception to this is that the orbitals tend to be more 
compact when <j > and that the orbitals tend to be 
more diffuse when a < than they are when cr = 
(no lattice relaxation). Tables 2 and 3 contain the 
numerical predictions of the point-ion model for U 
centers in KC1 and CaF-2 when the lattice relaxation is 
zero (o" = 0). These two tables show that the ordering 
of the singlet energy levels with increasing energy is 
l S(ls, Is), ^(ls, 2p), and l S(ls, 2s). They also show that 
the energies of the triplet states 3 S(ls, 25) and 3 P(ls, 
2p) lie between the ground singlet state energy and 
the first-excited-singlet state energy. The energies of 
these two low-lying triplet states lie close together and 
may be nearly degenerate as in the case for KC1. In 
addition the ordering of these two triplet energy levels 
depends upon the crystal and upon the lattice relaxa- 
tion. The results of [17] for the alkali halides predict 
the same ordering for the singlet energy levels. They 
also predict that the triplet 3 S(ls, 25) and 3 P(ls, 2p) 
energy levels lie between the l S(ls, Is) and l P(ls, 2p) 
energy levels. However, contrary to the present results, 
the results of [17] predict that the 3 S(ls, 2s) energy 
level lies below and is well spearated from the 3 P(ls, 
2p) energy level. These differences among the ordering 
of variational states and numerical triplet states may 
be explained by either or both of two possibilities: 
First, the authors of [17] have estimated the exchange 
terms in a different manner than the present author. 
And second, the numerical wave functions are more 
"flexible" than two-parameter-variational wave func- 
tions. One measure of the flexibility of a wave function 
is its second order-spatial derivative. The latter is 
proportional to the kinetic energy. Variational functions 
introduce artificial bounds on the kinetic energy and 
for some self-consistent potentials predict lower energy 
levels and different orderings of the states than 
numerical functions. 

Table 4 gives for the six crystals the predicted 
energy levels, f/-band transition energies, the [/-band 
oscillator strengths, and the lattice relaxation of the 
nearest-neighbor ions. The oscillator strengths / in 



TABLE 2. The five lowest-lying states for KC1 of the U center in the 
perfect point-ion lattice (cr = 0) with exchange energy and Coulomb 
correlation energy. For convenience, the states l S{\s, Is), : 'P(ls, 
2p), 3 S(Lj, 25), ! P(l5, 2p), and l S{ls, 2s) are labeled A, B, C, D, and 
E, respectively, in this table. The total energy of state X is Et (X; 
cr = 0), where X = A, B, C, D, or E. The transition energy from 
stated to state Y is E(X, Y). The quantities n and / are the prin- 
cipal and angular-momentum quantum numbers for the symmetry 
of the one-electron orbitals from which the {/-center configuration 
is made. The spatial extent quantities r„i(s) and r,(nl) are dimen- 
sionless and the energies are expressed in terms of atomic units. 
(1 a.u. of energy = 27.2 eV = 4.36 attojoules.) 



State 


A 


B 


C 


D 


E 


E T 


-1.015 


-0.953 


-0.957 


-0.878 


-0.856 


r«(D 


0.565 


.281 


.285 


.479 


.531 


r,„(3) 


.334 


.057 


.059 


.234 


.293 


r e (nl) 


.590 


.202 


.208 


.490 


.551 


AVr(l) 


.565 


.857 


.882 


.881 


.956 


iv ,.(3) 


.334 


.868 


1.005 


.929 


1.915 


re(n'l') 


.590 


1.012 


1.140 


1.054 


2.002 


E(A, D) 


.137 










E(D, E) 


.022 










E(B, C) 


.004 











Table 3. The five lowest-lying states for CaF-2 of the U center in 
the perfect point-ion lattice (cr = 0) with exchange energy and 
Coulomb correlation energy. The notation in this table is the 
same as the notation given in table 2. 



State 


A 


B 


C 


D 


E 


E T 


-1.204 


-1.072 


-1.047 


-0.999 


-0.960 


r„,0) 


0.626 


0.374 


0.372 


.559 


.553 


r„i(3) 


.435 


.123 


.120 


.340 


.330 


re(nl) 


.694 


.328 


.323 


.607 


.597 


rvr-(l) 


.626 


.940 


1.079 


.986 


1.516 


rvr(3) 


.435 


1.268 


2.180 


1.455 


1.608 


rAn'l') 


.694 


1.350 


2.021 


1.476 


4.012 


E(A, D) 


.205 










E(D, E) 


.039 










E(B, C) 


.025 











TABLE 4. Energy levels and (7-band transition energies for absorp- 
tion from the l S(ls, Is) state to the ! P(l5, 2p) state for the (7 center. 
The lattice relaxation and Coulomb correlation are included. The 
initial state is ^(ls, Is) and cr = en. The value of cr remains the 
same for the final state ' P(ls, 2p). For convenience, the states 
! S(b, Is) and ^(ls, 2/>) are labeled A and D, respectively, in this 
table. The total energy of the state XisE T {X; cr), where A' = ^4 or D. 
The (7-band energy is E(A, D). The experimental value for the 
(7-band energy is E(A, D; expt). The theoretical oscillator strength 
is/. The quantities cr, and /are dimensionless. All other quantities 
are in atomic units. (1 a.u. of energy = 27.2 eV = 4.36 attojoules.) 



Sub- 
stance 


cr, 


/ 


E T (A,<n) 


£/(/>; o-,) 


E(A,D) 


E(A,D, 

expt) 


NaCl 


0.041 


1.90 


-1.087 


-0.914 


0.173 


a 0.238 


KC1 


.060 


1.90 


-1.023 


-.870 


.153 


a .213 


CdF 2 


.050 


1.94 


-1.221 


-.987 


.234 




CaF 2 


.060 


1.94 


-1.221 


-.985 


.236 


b .281 


SrF 2 


.057 


1.94 


-1.177 


-.960 


.217 


b .259 


BaF 2 


.065 


1.94 


-1.138 


-.937 


.201 


b .221 



a J. W. Schulman and W. D. Compton, Color Centers in Physics 
(The Macmillan Company, New York, 1962). 

b J. H. Beaumont et al., J. of Physics C 3, L153 (1970). 
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table 4 agree reasonably well with other theoretical 
calculations [1, 16,21]. 

Referring to tables 4 and 5, we observe that including 
the lattice relaxation of the nearest-neighbor ions 
increases the predicted [/-band transition energies. 
Nevertheless, the point-ion model with lattice relaxation 
still yields values for the f/-band energies which are 
smaller than the observed values. Table 6 shows that 
these results are consistent with the other results 
predicted by point-ion models based upon variational 
wave functions. 

One possible explanation for the energy levels 
of the singlet states in the present work lying above 
the corresponding energy levels in variational theories 
is given in [21]. 

Tables 7 and 8 give the energy levels and transi- 
tion energies for absorption respectively from the x P(l5, 
2p ) state to the *S ( Is , 25 ) state and from the 3 P ( Is , 2p ) 
state to the 3 S(ls, 2s) state. 



Table 5. Energy levels and U-band transition energies for absorp- 
tion from the l S{ls, Is) state to the ^(ls, 2p) state for the U 
center with no ionic polarization. Because (7 = 0, the lattice 
relaxation is zero. The Coulomb correlation is included. The 
initial state is x 5(l5, Is) and C7=(7i = 0. For convenience, the 
states 1 S(ls, Is) and x P(l5, 2p) are labeled A and Z), respectively, 
in this table. The total energy of the state X is E T (X; cr), where 
X — A or D. The U-band energy is E(A, D). The experimental 
value for the c7-band energy is E(A, D; expt). The quantity cr, is 
dimensionless. All other quantities are in atomic units. (1 a.u. of 
energy = 27.2 eV = 4.36 attojoules.) 



Substance 


<Ti 


E T (A;(Ti) 


E T (D;<Ti) 


E{A,D) 


£ (/*, D; expt) 


NaCl 


0.0 


- 1.082 


-0.921 


0.161 


a 0.238 


KC1 


.0 


-1.015 


-.878 


.137 


a 0.213 


CdF 2 


.0 


-1.211 


- 1.003 


.208 




CaF 2 


.0 


- 1.204 


- 0.999 


.205 


b 0.281 


SrF 2 


.0 


- 1.162 


- 0.972 


.190 


b 0.259 


BaF 2 


.0 


-1.117 


- 0.943 


.174 


b 0.221 



Table 6. Comparisons among various theoretical values and experimental values for the U-band transition energies. This is the transition 
from the l S(ls, Is) state to the 1 P(ls, 2p) state. All energies are in atomic units (1 a.u. of energy = 27.2 eV = 4.36 attojoules). 

















Present 


Substance 


Experiment 


Ref. a 


Ref. b 


Ref. c 


Ref. d 


Ref. e 


work 
(Table 4) 


NaCl 


f 0.238 






0.203 


0.202 




0.173 


KC1 


f .213 


0.205 


0.219 


.178 


.178 




153 


CdF 2 














.234 


CaF 2 


g .281 










0.251 


.236 


SrF 2 


*.259 










.234 


.217 


BaF 2 


K .221 










.217 


.201 



a R. F. Wood and U. Opik, Phys. Rev. 162, 736 (1967); first model. 

b R. F. Wood and U. Opik, Phys. Rev. 162, 736 (1967); second model. 

c H. N. Spector et al., J. Chem. Phys. 46, 2676 (1967). 

d B. S. Gourary, Phys. Rev. 112, 337 (1958). 

e R. S. Singh et al., J. Chem. Phys. 52, 2341 (1970). 

f J. W. Schulman and W. D. Compton, Color Centers in Physics (The Macmillan Company, New York, 

B J. H. Beaumont et al., J. of Physics C 3, L153 (1970). 



1962). 



a J. W. Schulman and W. D. Compton, Color Centers in Physics 
(The Macmillan Company, New York, 1962). 

b J. H. Beaumont et al., J. of Physics C 3, L153 (1970). 



Table 7. Energy levels and second singlet transition energies for 
absorption from the ! P(ls,2p) state to the l S(ls, 2s) state for the 
U center. The lattice relaxation and Coulomb correlation are 
included. The initial state is l P(ls, 2p) and cr = cr/. The value of 
cr remains the same for the final state *S(ls, 2s). For conven- 
ience, the states ^(ls, 2p) and 'Sds, 2s) are labeled D and E, 
respectively. The total energy of the state X is Et(X; cr). where 
X=D or E. The second singlet transition energy is E(D, E). The 
quantity cr, is dimensionless. All other quantities are in atomic 
units. (1 a.u. of energy = 27.2 eV = 4.36 attojoules.) 



Substance 


cr, 


E T (D; en) 


E T (E; cr,-) 


E(D.E) 


NaCl 


-0.011 


- 0.921 


- 0.894 


0.028 


KC1 


.000 


-.878 


-0.856 


0.022 


CdF 2 


-.082 


-1.010 


- 1.010 


0.000 


CaF 2 


.000 


- 0.999 


- 0.960 


0.039 


SrF 2 


.001 


-.972 


- 0.935 


0.037 


BaF 2 


.024 


.945 


- 0.904 


0.042 



Table 8. Energy levels and first triplet transition energies for 
transitions from the 3 P(Ls, 2p) state to the 3 S(ls, 2s) state for 
the U center. The lattice relaxation and exchange energy are in- 
cluded. The initial state is 3 P(ls, 2p) and cr = cr,. The value of cr 
remains the same for the final state 3 S(ls, 2s). These states are 
labeled B and C, respectively, for convenience. The total energy 
of the stated is Et(X, cr), where ^ = # orC The first triplet transi- 
tion energy is E(B, C). The quantity cr, is dimensionless. All 
other quantities are in atomic units. (1 a.u. of energy = 27.2 
eV = 4.36 attojoules.) 



Substance 


cr, 


E T (B; en) 


E T {C; en) 


E{BX) 


NaCl 


- 0.006 


- 0.996 


- 0.992 


0.004 


KC1 


0.004 


- 0.953 


- 0.957 


- 0.004 


CdF 2 


- 0.014 


- 1.077 


- 1.057 


0.021 


CaF 2 


0.009 


- 1.072 


- 1.044 


0.029 


SrF 2 


0.006 


- 1.046 


- 1.026 


0.020 


BaF 2 


0.028 


- 1.020 


- 1.004 


0.017 
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A Mollwo-Ivey relation [24J. 

E = aa~ r \ 

with n ~ 1.1 relates the observed energies E at 
the peak of the (/-band in many alkali halides to 
the lattice spacing a. However, the observed energies 
of the [/-band peaks in the alkaline earth fluorides do 
not fit a Mollwo-Ivey empirically-derived ralation [20]. 
A similar situation occurs for the F-bamd peaks. A 
Mollwo-Ivey relation with n ~ 1.84 predicts quite 
well the locations of the F-band peaks in many NaCl- 
type alkali halides. But, such a relation does not pre- 
dict the observed F-band peaks in the alkaline earth 
fluorides [20]. 

Point-ion models for the (/-band which do not include 
lattice relaxation yield (/-band energies which follow 
Mollwo-Ivey relations. The variational-point-ion 
theory of [17] predicts (/-band energies for NaCl- 
type crystals which obey a Mollwo-Ivey relation with 
n ~ 1.14. The similar variational theory of [1] gives 
(/-band energies for the alkaline earth fluorides which 
follow reasonably well a Mollwo-Ivey relation with 
n ~ 1.143. The present point-ion model in which the 
lattice relaxation and the coulomb correlation term 
W(r) are zero predicts that the Mollwo-Ivey exponents 
are 1.22 for (/-bands in NaCl and KC1 and 1.31 for (/- 
bands in the fluorides CdF 2 , CaF 2 , SrF2, and BaF2. 
When the Coulomb correlation term W(r) in a point- 
ion model with no lattice relaxation is included, the 
Mollwo-Ivey exponents become, respectively, 1.46 
and 1.33. The predicted Mollwo-Ivey exponents for 
the alkali halides which are obtained by variational 
methods are closer to the experimental values than 
those exponents which are obtained by the present 
numerical methods. However, when the lattice 
relaxation is included, then the predicted (/-band 
energies for both the alkali halides and the fluorides 
do not follow Mollwo-Ivey relations. 

The above results show that solving the HFS 
equations for the point-ion model of the U center by 
numerical techniques leads to a greater disagreement 
between theory and experiment for the (/-band ener- 
gies than does solving the point-ion model by variational 
techniques. Also, the Mollwo-Ivey exponent obtained 
from the numerical HFS wave functions differs more 
from the observed exponent for (/-bands in the alkali 
halides than does the exponent obtained from the vari- 
ational wave functions. These greater disagreements 
between the predictions of the numerical HFS wave 
functions and predictions of the variational wave func- 
tions for the point-ion model suggest that the past suc- 
cesses of the point-ion model have been fortuitous. 
More generally, the disagreement among the observed 
and computed values reported in this paper lead one 
to conclude that the point-ion models do not represent 
satisfactorily the properties of the (/-center. The pres- 
ent work also emphasizes a need for models containing 
more detailed treatment of ionic polarization, electronic 
polarization, and the finite extent of the neighboring 
ions. Examples of two such models for the alkali halides 
appear in [18]. 



Whenever both of the (/-center electrons have com- 
pact wave functions, the U center (proton, anion 
vacancy and two electrons) is essentially neutral. 
In this case, it would produce only small, long-range 
polarizations, if any. The existence of these conditions 
then would justify treating all ions as point charges. 
Hence, the point-ion model is a reasonable approxi- 
mation for the case in which most of the charge associ- 
ated with the (/-center electrons is localized about the 
proton-anion vacancy site and remains outside the 
neighboring ionic cores. The defect electrons then 
would experience a potential which scales with the 
lattice constant and the predicted U bands would vary 
with the lattice constant, to a good approximation, 
according to a Mollow-Ivey relation. 

The observed U bands in the alkali halides satisfy 
a Mollwo-Ivey relation; those in the alkaline earth 
fluorides do not. The failure of the observed U bands 
in the alkaline earth fluorides to obey a Mollwo-Ivey 
relation suggests the speculation that extended ion 
effects are more important in the alkaline earth 
fluorides than they are in the alkali halides. Extended 
ion effects include exchange and overlap among the 
(/-center electrons and the core electrons of neigh- 
boring ions. These effects probably do not depend 
smoothly upon the lattice constant. It is interesting to 
notice in table 6 that most of the predicted U bands 
in the alkaline earth fluorides are closer to the observed 
values than the predicted U bands in the alkali halides. 
Nevertheless, using a point-ion model to describe the 
U center in the alkaline earth fluorides is more ques- 
tionable than using a point-ion model to describe the 
U center in the alkali halides because the former class 
of U centers fails to follow a Mollwo-Ivey relation. 
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4. Appendix. Assumptions Contained in the 
HFS Equations 

Because the Schrodinger equation for the two-elec- 
tron U center, 

H T (x, y; o-) ¥„(*, y)=EVr,(x, y), 

cannot be solved exactly, we use the self-consistent 
field method [Hartree-Fock (HF)] to calculate the wave 
functions of the U center. The wave functions for the 
self-consistent field approximation are antisymmetrized 
sums of products of one-electron functions Ui(r) and 
u 2 (r): 

^(x, r )-^(*, y;HF) 

= 2~ 1 > 2 [u l (x)u 2 (y)-u,(y)u 2 (x)]. 
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The functions u have a restricted-central-field rep- 
resentation given by 

Unlms=Rnl(r)yim(0, (f))a s , 

where R n i(r) is the radial function, Yi m is a, spherical 
harmonic function, and a s is the spin function. The 
principal, orbital, magnetic, and spin quantum num- 
bers are, respectively, n, /, m 9 and s. The radial func- 
tions Pni(r), P n i(r) = rR n i(r), are normalized to the 
crystal volume f P n i 2 (r) dr=l. The spherically aver- 
aged total electronic charge density for both spins in 
terms of the radial functions becomes 

Pav(0=-e[o-(r)/477r 2 ], 

where the spherical density cr(r) is 

o-(r) = ^ (DniPl,{r), 
nl 

and the occupation number of the spatial orbital for 
both spins is a> n i. The summation S«/co„/ equals 2 for 
the U center. 

In terms of the above assumptions, the Coulomb 
potential V c (r) for the HFS equations then has the 
form: 

V c (r)=- I' <r(t)dt + e 2 f*^-dt. 

r JO Jr t 

Slater [23] suggests that the essential features of the 
HF method are retained when a common exchange 
potential replaces the exchange potentials for each of 
the orbitals U\ and u 2 . Applying his suggestion to the 
triplet states of the U center, we assume that the com- 
mon exchange potential V e (r) for the U center at the 
point r is equal to the exchange potential which a 
free-electron gas would have if its total electron charge 
density for both spins were equal to that of the non- 
uniform system (the two [/-center defect electrons): 
namely, 



V e (r)=Se 2 [(3/8tt) e-'\p„ v (r)\] 



1/3 



Let us define V s (r) = Vc(r) + V e (r) for future use. 
The self-Coulomb energy part of the total electronic 
Coulomb energy cancels exactly in the conventional 
HF method a corresponding energy in the total ex- 
change energy. This cancellation does not occur in the 
HFS equations at large r. We alter the sum V 8 (r) so 
that it has the correct asymptotic behavior at large r. 
We define V(r) = V c (r) + V e (r) for r < r and V(r) 
= e 2 {XniO) n i — l)/r for r^ r () . The radius r is that value 
of r at which 



V s (ro) = e 2 (^com=l)lro. 

\ nl I 



The HF method also does not include the spatial 
correlation in the motion of the two defect electrons 
produced by their instantaneous Coulomb repulsion 



W{r)=-e 2 



0.288 



5.1 a + r s (r)' 



0^W(r 8 )IV exch (r 8 ) 



0.314r,/(5.1ao+r s )^0.314r s 



7/ r (x, x). The approximate HF wave function is the 
source of this Coulomb correlation problem because it 
does not depend upon the distance between the two 
electrons |x — y|. We define the correlation energy by 
the relation 

E e - <¥,(*, y)|»r|¥,(*,y)> 

-<**(*, r; HF) \H T \Vr,(x, y; HF)). 

Let us assume also that a correlation operator € c ( |x 
— y | ) exists such that its expectation value in the 
approximate representation 4s(x, y; HF) is the 
correlation energy 

J? c = <¥„(*, y; HF\€ C ( |x-y|) |<M*. y; HF)). 

We may introduce, then, the correlation potential 
W c {r) by the following operation: 

W c {r) = JdHVvix, r; HF)e c { \ x-r| )¥»(*, r; HF). 

Mitler [22] adapts to atomic helium the free-electron 
gas approximation for the correlation energy. He 
introduces in the HF equations an additional central 
"correlation" potential W(r) to which pairs of electrons 
with opposite spin are subject. The approximate cor- 
relation potential W (r) has the form [22] 



where a is the Bohr radius and the local density of 
electrons is 

e" 1 |p('-)l=[(477/3)r 2 s (r)]-'. 

We expect that the expression for W(r) is a good 
estimate for the correlation potential W c (r). Mitler 
[22] applies the expression for W{r) to the ground state 
for helium and obtaines by perturbation theory a 
ground state energy which agrees to within J percent 
of the experimental value. Our prescription is, then, 
to replace V c (r) in the singlet HFS variational equa- 
tions with V c (r) + W(r). Because the inequality 



is obtained, we do not include the Coulomb correlation 
potential in the triplet HFS variational equations. 
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